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ABSTRACT 

Observations show that star formation is an inefficient and slow process. This result can be at- 
tributed to the injection of energy and momentum by stars that prevents free-fall collapse of molecular 
clouds. The mechanism of this stellar feedback is debated theoretically: possible sources of pressure 
include the classical warm HII gas, the hot gas generated by shock- heating from stellar winds and su- 
pernovae, direct radiation of stars, and the dust-processed radiation field trapped inside the HII shell. 
In this paper, we measure observationally the pressures associated with each component listed above 
across the giant HII region 30 Doradus in the Large Magellanic Cloud. We exploit high-resolution, 
multi-wavelengh images (radio, infrared, optical, and X-ray) to map these pressures as a function of 
position. We find that radiation pressure dominates within 75 pc of the central star cluster, R136, 
while the HII gas pressure dominates at larger radii. By contrast, the dust-processed radiation pres- 
sure and hot gas pressure are generally weak and not dynamically important, although the hot gas 
pressure may have played a more significant role at early times. Based on the low X-ray gas pressures, 
we demonstrate that the hot gas is only partially confined and must be leaking out the HII shell. 
Additionally, we consider the implications of a dominant radiation pressure on the early dynamics of 
30 Doradus. 

Subject headings: galaxies: star clusters — HII regions — ISM: individual (30 Doradus) — stars: 
formation 



1. INTRODUCTION 

Molecular clouds contain the coolest and densest gas 
in the Universe, and thus they are the sites where stars 
form. The physical properties of these clouds set the ini- 
tial conditions for protostellar collaps e and may define 
the stellar initial mass function (IM F) iMot te et al .11 1998c 
iTesti fe Sargentl 11998c lOnishi et all [2002[. The massive 
stars formed there eventually end in supernova explo- 
sions, injecting mechanical energy and chemically enrich- 
ing the interstellar medium (ISM). Therefore, molecular 
clouds shape the entire stellar life cycle, and an under- 
standing of their properties and dynamics is key to probe 
galactic evolution. 

Observational evidence shows that star formation is 
an inefficient and slow process. Only 5-10% of avail- 
able molecular clou d mass is converted into stars over 
the cloud lifetimeS (jWilliams fe McKed fl997h and only 
^2% of the gas is converted to stars in one free- 
fall time across several orders of magnitude in den - 
sity (jZuckerman fc Evans! fl974c IKrumholz fe Tanll2007D . 
This inefficiency can be attributed to the internal pro- 
cesses of HI I regions that disrupt their host mole cular 
clouds (e.g., iMatznerl l200l IKrumholz et al.ll2QQ6h . but 
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the mode of this stellar feedback remains uncertain. 

Broadly, there are several possible sources of inter- 
nal energy and momentum that may drive the dy- 
nami cs of HII regions: the direct radiation from star s 
(e.g.. Uijina fc Adamsl[T996c IKrumholz fc Matznerll2009[ ). 
the dust-p rocessed infrared radia t ion trapped inside an 
HII shell (iThompson etlll [2005c iMurrav et all l2010bc 
lAndrews fc T hompson 2Q TTh . the warm gas ionized b y 
massive stars (e.g.. IWhitworthl 11979b iDale et al.ll2QQ5h . 
the hot gas shock heated by stellar winds and super- 
novae (e.g., lYorke et al.l 1 1989c lHarper-Clark fc Murray! 
2009) , and protostellar ou tflows/ jets (e.g.. iQuillen et al.l 
2005; iCunnin gham e faH 120061; iLi fc Nakamural [2005 
iNakamura fc Lill2008c (Wang et al.ll2010D . Each of these 
mechanisms has been considered individually in the liter- 
ature, but no observational analyses have ever compared 
the relative contribution of all these components within 
HII regions. 

In this paper, we investigate the role of the stellar feed- 
back mechanisms listed above in the giant HII region 30 
Doradus in the nearby Large Magellanic Cloud (LMC). 
Several properties of the LMC make it a favorable target: 
the LMC's proximity (~50 kpc) ensures individual point 
sources can be resolved while maintaining the capability 
of mapping the diffuse emission at sub-pc scales. Addi- 
tionally, the LMC has a face-on orientation and a low 
column density (a few xlO 21 cm -2 ) that limits line-of- 
sight confusion. Given these advantages, the LMC (and 
thus 30 Doradus) has been surveyed at many wavelengths 
at high spatial resolution, and we can exploit these data 
to compare observationally all the feedback mechanisms 
and how they vary with position across 30 Doradus. 

The text is structured as follows: §1.1 gives relevant 
background on the source, 30 Doradus, and describes 
why this source is a good "test case" for our analyses. 
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In §2, we present the multiwavelength data utilized in 
our work to assess the dynamical role of all the possible 
stellar feedback mechanisms. §3 outlines how we utilize 
these images to calculate the pressures associated with 
each feedback component across 30 Doradus. §4 gives 
the results from our analyses, and §5 discusses the im- 
plications of our findings, including evidence of X-ray 
gas leakage from the HII region (§5.1) and the role of 
radiation pressure in HII region dynamics (§5.3). Addi- 
tionally, we articulate the different ways one can define 
radiation pressure, and how these definitions can lead 
to divergent results in §5.2. Finally, we summarize and 
conclude our analysis in §6. 

1.1. Background on 30 Doradus 

30 Doradus is the most massive and largest HII region 
in the Local Group. The primary star cluster powering 
30 D oradus is NGC 2070, with 2400 OB stars (jParkerl 
fl993h . an ionizing photon luminosity S' = 4.5xl0 51 pho- 
tons s _1 (Walborn 1991), and a bolometric luminosity of 
7.8xlO 7 L (IMalumuth fe Head[l99l. The IMF of NGC 



1(f 



2070 has masses up to 120 M (jMassev fc Hunterll 19981 ). 
and the stellar population may be the result of se veral 
epochs of star formation (Walborn & Blades 1997). At 
the core of NGC 2070 is R136, the densest concentration 
of very massive stars known, with a central density of 5.5 
xlO 4 M pc~ 3 (jHunter et all 11995ft : R136 hosts at least 
39 03-type stars and at least 10 WR stars in its core 
(~2.5 pc in diameter; iMassev fc Hunted [199 8). 

To provide context for how 30 Doradus compares to 
other local HII regions, Figure [1] plots Ha luminosity 
versus HII region radius for ~22,000 HII regions in 70 
nearby (distances <30 Mpc) galaxies (see references in 
figure caption). Morphologically, this galaxy sample is 
comprised of 13 irregulars/dwarf irregulars and 57 spi- 
rals. The black star near the top right denotes 30 Do- 
radus. It is the brighest in Ha of the 613 HII regions in 
the irregulars by nearly an order of magnitude. Relative 
to the HII regions in spirals (including M33), 30 Doradus 
has a greater Ha luminosity than ~99% of that sample. 

The nebula that surrounds the central star cluster has 
a complex morphology across the electromagnetic spec- 
trum. Figure [2] shows a three-color image of 30 Doradus, 
with the Spitzer Space Telescope 8-/im IRAC band in 
red, Ha in green, and soft X-rays (0.5-2.0 keV) in blue 
(details of these data are given in §2). Large and small- 
scale structures are evident, from thin ionized gas and 
dust filaments of arcsecond widths to cavities a few ar- 
cminutes across filled with hot X-ray gas. The warm 
ionized gas has several shell-like structures, and many of 
these are expanding w ith high velocities (~ 100-300 km 
s -1 : lChu fc Kennicuttlll994f ) . suggesting that past super- 
nova explosions have occurred in the region. In addition 
to a large ionized gas mass (~ 8 x 1O 5 M ; iKennicuttl 
11984ft. the 30 Doradus nebula also has - 1O 6 M of CO 
(j Johansson et all 11998( 1. The CO(1-0) maps of 30 Do- 
radus have revealed 33 molecular cloud complexes in the 
HII region, and in particular, two elongated clouds of 
CO mass ~ 4 x 1O 5 M that form a "ridge" West and 
North of R136 (see the CO contours in Figure [2]). Esti- 
mate s of the radi us Rmi of the nebula range from ~110 
pc (jBran dl 2005; using a rev ised value of D = 50 kpc) 
to ~185 pc (jKennicuttl Il984f ) . The nearly factor of two 
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Fig. 1. — Ha luminosity versus HII region radius for ~22,000 
HII regions in 70 nearby (distances < 30 Mpc) galaxies. The black 
star at the top right denotes 30 Doradus. It is the brightest 
HII region in the irregular galaxies by nearly an order of mag- 
nitude, and it is more luminous than ^99% of the HII regions in 
the spiral galaxies. We note that none of the data is corrected 
for reddening. The plotted dat a were compiled from the follow- 
ing refe rences: LMC and SM C: [Kennicutt & Hodge (1986) ; Sex - 
tans A: Hodg elFaLl (1199431): NG C 6822: [Hodge et al. (1989b) 
Holmberg II: IHodge et alJ (11994b): GR8: Hod ge et al.l fl989a) 
DDO 47, Le o A, S extans B, DDO 167, DDO 168 , DDO 187 
I Strobel et all (fl99ll) : DDO 53 Ejtr obel et al. (1990); 56 spirals 
Knapen et al. (2003), Bradle y et al.l (|2006l ): M33: IWvder et all 
(1997), Hodge et al. (1993). We utilized the distances from 
Kennicutt et al. (2008) to convert from Ha flux to luminosity. 

uncertainty in i?Hii arises from the complex shape that 
precludes accurate determination of the radius. In this 
paper, we will assume .Rhii = 150 pc. 

The properties of 30 Doradus described above demon- 
strate why this HII region is an ideal candidate for assess- 
ing the feedback mechanisms of massive stars. The shear 
number and energetic output of the OB stars facilitates 
detailed study of the effects of radiation, winds, super- 
novae, ionization fronts, etc. Additionally, the proximity 
of 30 Doradus enables a resolved view of the processes 
and dynamics associated with starburst activity that 
was c ommon in the early Universe (e.g., Me urer et al.l 
119971 : iShaplev et all l2QQ3h . Indeed, the relatively in- 
stantaneous formation of the concentrated ma ssive stars 
in R1 36 make 30 Doradus a "mini-starburst" (jLeithererl 
fl997h . 

2. DATA 

We analyzed images of 30 Doradus at several wave- 
lengths. A brief description of these data is given below. 

2.1. Optical 

We compiled optical photometric data on 30 Do- 
radus from three separate observational programs. For 
the central 35" x 35"around R136 (with right ascen- 
sion a = 05 h 38 m 45.5 s and declina tion -69°06 / 02.7 // ) 
we ut ilize the photometric results of IMalumuth fc Heap 
(|1994f ) from Hubble Space Telescope Planetary Camera 
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Fig. 2. Three-color image of 30 Dor: MIPS 8/xm (red), Ha 
(green), and 0.5-8 keV X-rays (blue). White contours show the 
12 CO(1-0) emission (Johansson et al. 1998) in the region. Both 
large- and small-scale structures are evident. North is up, East is 
left. 



observations. These authors identified over 800 stars 
within this area and obtained a bolometric luminosity 
^boi = 7.8 x 1O 7 L for their sources. 

At larger distances from R136 out to a f ew arcminutes , 
we employ the UBV photometric data of Parker (1993). 
The optical images of 30 Doradus from Parker (1993) 
were obtained at the 0.9-m telescope at Cerro Tololo 
Inter- American Observatory (CTIO), with a field of view 
of 2.6' x 4.1' and O^'pixel - 1 . We followed the analyses 
of iParker fc Garmanvl (|1993f ) to convert their measured 
apparent UBV magnitudes to absolute bolometric mag- 
nitudes. 

For the area out side the field of [ Parker (1993), we use 
the UBV data of lSelman fc Melnickl (|2005j). These ob- 
servations were taken with the Wide Field Imager on the 
MPG/ESO 2.2-m telescope at La Silla, out to half a de- 
gree away from R136 with 0.238"pixel -1 . Thus, the three 
datasets combined provide full coverage of 30 Doradus in 
the U, B, and V bands. 

To illustrate the HII region structure, we show the 
Ha emission of 30 Doradus in Figure [2j This narrow- 
band image (at 6563A, with 30A full-width half max 
was taken with the University of Michigan/ CTIO 61- 
cm Curtis Schmidt Telescope at CTIO as part of 
the Magellanic Cloud Emissi on Line Survey (MCELS: 
ISmith fc MCELS Teaml [1998). The total integration 
time was 600 s, and the reduced image has a resolution 
of 2"pixer 1 . 

2.2. Infrared 

Infrared images of 30 Doradus were obtained through 
the Spitzer Space Telescope Legacy project Surveying 
the A gents of Galaxy Evolution (SAGE: Mei xner et al.l 
2006) of the LMC. The survey covered an area of —7 
x 7 degree s of the LMC wit h the Infrared Array Cam- 
era (IRAC; iFazio et al.l 12004ft and the Mu ltiband Imag- 
ing Photometer (MIPS; iRieke et a"T1l2004l ). Images were 
taken in all bands of IRAC (3.6, 4.5, 5.8, and 7.9 /mi) 
and of MIPS (24, 70, and 160 /im) at two epochs in 2005. 
For our analyses, we used the combined mosaics of both 



epochs with l^'pixel 1 in the 3.6 and 7.9 /im IRAC im- 
ages and 2.49 // pixer 1 and 4.8 // pixels~ 1 in the MIPS 24 
/im and 70 /im, respectively. 

2.3. Radio 

30 Doradus was observed with the Australian Tele- 
scope Compact Array (ATCA) as part of a 4 . 8- GH z 
and 8.64-GHz survey of the LMC (jDickel et al.ll2005l ). 
This program used two array configurations that pro- 
vided 19 antenna spacings, and these ATCA observa- 
tio ns were combined wi th the Parkes 64-m telescope data 
of lHavnes et al.l (|1991[ ) to account for extended struc- 
ture missed by the inter ferometric observations. For our 
analyses, we utilized the resulting ATCA+Parkes 8.64 
GHz (3.5-cm) image of 30 Doradus, which had a Gaus- 
sian beam of FWHM 22" and an average rms noise level 
of 0.5 mJy beam -1 . We note that higher-resolution 
ATCA observa ti ons o f 30 Doradus have been taken by 
lLazendic et al.l (|2003[ ). but we have opte d to use the 
ATCA+Parkes image of lDickel et~aT1 (|2005f ) as the latter 
is more sensitive to the low surface-brightness outskirts 
of 30 Doradus. 

2.4. X-ray 

30 Doradus was observed using the Chandra Advanced 
CCD Imaging Spectrometer (ACIS) in 2006 January for 
^94 ks total (ObsIDs 5906 [13 ks], 7263 [43 ks], 7264 
[38 ks]; PI: L. Townsley) in the Timed- Exposure VFaint 
Mode. The spatial resolution of the Chandra ACIS im- 
ages is 0.492" pixel -1 . Data reduction and analysis was 
performed using the Chandra Interactive Analysis of Ob- 
servations (ciao) Version 4.1. We followed the CIAO data 
preparation thread to reprocess the Level 2 X-ray data 
and merge the three observations together. Figure [3] 
shows the resulting soft X-ray band (0.5-2.0 keV) im- 
age following these analyses. Seventy- four point sources 
were identified in the reprocessed images using the CIAO 
command wavdet ect (a source detectio n algorithm using 
wavelet analysis; iFreeman et al.l 12002); we excluded the 
identified point sources in our spectral analyses. 

To produce a global X-ray spectrum of 30 Doradus, 
we extracted Chandra spectra using the CIAO command 
specextract. Background spectra were also produced 
from a circular region of radius ^15" that is ~2' East 
of 30 Doradus, and these were subtracted from the 
source spectra. Additionally, we removed the counts 
of the 74 point sources identified above. The result- 
ing spectra were modeled simultaneously as an absorbed, 
variable- abundance plasma in collisional ionization equi- 
librium (XSPEC model components phabs and vmekal) 
in XSPEC Version 12.4.0. Figure H gives the spectra 
with the best-fit model (with y 2 = 619 with 396 degrees 
of freedom [d.o.f.]) overplotted. We found a best-fit ab- 
sorbing column density of TVh = 1-51q"J x 10 21 cm -2 and 
an X-ray gas temperature of kTx = 0.641q q2 keV. The 
absorption-corrected soft-band (0.5-2.0 keV) luminosity 
of the diffuse emission is Lx = 4.5 x 10 36 erg s -1 . 

Previou s Chandra X-ray analysi s of 30 Doradus was re- 
ported bv lTownslev et al] (2006a. b) for a different set of 
observations (ObsIDs 22 and 62520) totalling -24 ks. By 
fitting the X-ray spectra of many diffuse regions across 
30 Doradus, they found best-fit absorbing columns of 
N u = 1 - 6 x 10 21 cm -2 , temperatures of kT x ~ 0.3-0.9 
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Fig. 3. Chandra ACIS X-ray soft band (0.5-2.0 keV) image of 
30 Doradus. The image was binned by a factor of four so structures 
are visually apparent. The cyan circle w 2' East of 30 Doradus is 
the area where background spectra were extracted. North is up, 
East is left. 
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Fig. 4. — The global X-ray spectra from the three ACIS observa- 
tions of 30 Doradus (ObsID 5906 in black, ObsID 7263 in red, and 
ObsID 7264 in green), with the best-fit models overplotted. The 
inset is the 68%, 90%, and 99% confidence contours for the column 
density Nn and the temperature kTx- The bottom panel gives the 
residuals between the data and the model in terms of x 2 . We find 
a best-fit N u = 1.5±° | x 10 21 cm" 2 and kT x = 0.64t°;^ keV. 

keV, and absorption-corrected luminosities (0.5-2.0 keV) 
of logLx = 34.2-37.0 erg s _1 . Thus, our values are fairly 
consistent with those of Townsley et al. 

3. METHODOLOGY 

To assess how feedback varies spatially across 30 Do- 
radus, we separate the source into 441 regions (see Fig- 
ure [5]). The area of the individual regions was selected to 
ensure sufficient signal-to- noise across the analyzed wave- 
bands; we chose the width of the regions (35"« 8 pc on 
a side, at a distance D = 50 k pc to the LMC) to matc h 
the HST PC image of R136 (Ma lumuth fc Heapl fl994). 
so that we could use their Lboi value and not have to 
resolve the individual point sources in the crowded R136 




Fig. 5. Ha image (from MCELS; Smith et al. 1998) with the 
441 regions we analyzed overplotted. Red squares denote those 
included in the hot gas leakage analysis of §5.1. The red X marks 
the center of R136. 

cluster. The number and position of our 441 regions was 
determined by the field-of-view and orientation of the 3- 
cm radio and Chandra X-ray data. Figure [5] shows the 
H-a image with all the resulting regions overplotted. 

To ascertain the dynamical importance of the feedback 
processes, we compute the pressures for each region us- 
ing the methods and relations described below. Since 
protostellar outflows a re only import ant dynamically in 
low-mass star clusters (jMatznerll2007l ). we do not expect 
them to play a big role in 30 Doradus, and we will not 
consider them in the rest of the text. 

3.1. Direct Radiation Pressure 

The light output by stars produces a direct radiation 
pressure that is associated with the photons' energy and 
momentum. The resulting radiation pressure Pdi r at 
some position within the HII region is related to the bolo- 
metric luminosity of each star Pboi and the distance r the 
light traveled to reach that point: 

where the summation is over all the stars in the field. 
In §5.2, we describe an alternative definition of radiation 
pressure and compare the results from each case. 

The above relation assumes that the stellar radiation 
are not attenuated by dust. In § 13. 2\ we calculate sepa- 
rately the radiation pressure associated with the light ab- 
sorbed by dust using Spitzer IR photometry. Given that 
our results show that Pdir ^> Pir generally (see §|4j), the 
assumption that the emitted Lboi is unattenuated seems 
reasonable. 

In order to obtain the bolometric luminosities of the 
massive stars in 30 Doradus, we utilize the UBV photo- 
metric data described in §2.1. To simplify the calcula- 
tion, we assume the bolomet r ic lu minosity of R136 ob- 
tained by iMalumuth fc Heapl (j!994f ) originates from the 
point in the middle of the central region marked with the 
red X in Fig. [5] For t he stars l ocate d outside R136 within 
a few arcminutes, the lParkerl 0993) catalog only includes 
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the apparent UBV magnitudes and colors. Therefore, 
we fo llow the procedure outlined by Parker & Gar manvl 



(|1993| ) to obtain absolute bolometric magnitudes of the 
126 4 stars in the Parker (1993) catalog that are not in 
the Selman & Melnick (2005) sample. For the 7697 stars 
in the iSelman fc Melnickl (120051) catalog that lie outside 
the field of the lPar ker (1993) data, we use their published 
values for the stars' absolute bolometric magnitudes. 

Thus, in total, we calculate the bolometric luminosi- 
ties Lboi of the R136 cluster and 8961 other stars in 30 
Doradus. For each of the 441 regions, we sum these 8962 
terms in Equation [TJ where r corresponds to the pro- 
jected distance from the 8962 stars' positions to the re- 
gion center. In this manner, we compute the radiation 
presure 'felt' by the 441 regions from all of the starlight 
in 30 Doradus. 

Since the stars are viewed in projection, the actual 
distance r to a star from the R136 center is observed 
as a projected distance ip. Therefore, we calculate the 
direct radiation pressure for two scenarios: one case as- 
suming the stars lie in the same plane (Pdir) and another 
case where we attempt to "deproject" the stars positions 
(^dir,3D)- Appendix A outlines the procedure we utilize 
to obtain the deprojected bolometric luminosity of the 
stars as a function of r and compares Pdi r with Pdir, 3D- 
We find that Pdir, 3D is 10-60% less than Pdi r at radii 
< 20 pc from 30 Doradus, and the fractional difference 
between Pdir, 3D and Pdi r at larger radii is much less (0.1- 
3.0%). As these differences do not affect the conclusions 
of this paper, we will only consider Pdi r for the rest of 
our analyses. 

3.2. Dust- Processed Radiation Pressure 

The stars' radiation will be processed by the nearby 
dust in the region, and an associated pressure is ex- 
erted by the resulting infrared radiation field, Pir. This 
pressure component could become dynamically impor- 
tant if the expanding HII shell is optically thick to the 
IR lig ht, effectively trapping the r adiation inside the HII 
shell (jKrumholz fc Matzn er 2009). The pressure of the 
dust-processed radiation field Pir can be determined by 
the energy density of the radiation field absorbed by the 
dust, u v (i.e., we assume steady state), 



Pir 



1 



(2) 



To find u v in each of our 441 regions, we measure 
their flux densities F v in the IRAC and MIPS images 
and c ompare them to th e predictions of the dust mod- 
els of iDraine fcTil 120071 (DL07 hereafter). The DL07 
models show how the IR emission spectral-energy dis- 
tribution varies depending on the dust content and the 
intensity of radiation heating the dust. DL07 assume 
a mixture of carbonaceous grains and amorphous sili- 
cate grains that have a size distribution that reproduces 
the wavel ength-dependent extinction in the local Milky 
Way (see IWeingartner fc Drain e 2001). One component 
of this dust mixture is polycyclic aromatic hydrocarbons 
(PAHs), small-sized carbonaceous grains that produce 
strong emission features at ~3-19 jam observed in many 
galaxies. 

Since the infrared emission from dust is relatively in- 
sensitive to the spectrum of the incident photons with 
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Fig. 6. — Measured IR flux ratios for the 441 regions (black 
diamonds) and the predicted flux ratios for different PAH mass 
fractions gpAH and scaling U of the energy density of the dust- 
processed radiation field (Equation [3} from DL07. We interpolate 
the grid of predicted flux ratios to obtain (/pah and U for each 
region. 

hv < 13.6 eV, DL07 adopts the spectrum of the local- 
neighborhood ISM. Then, u v is given by 



Uu 



IRSF 



(3) 



where U is a dimensionless scale factor and u][, RSF is the 
energy density of the hv < 13. 6 eV photons in the local 
ISM, 8.65 x 1(T 13 erg cm" 3 (jMathis et al.l[l983l ). We 
assume that each region is exposed to a single radiation 
field because the starlight heating the dust comes largely 
from NGC 2070. In DL07 parameters, this case corre- 
sponds to U m { n = /7 max and 7=0, where (1 — 7) is the 
fraction of the dust mass exposed to the radiation. 

For our analyses, we measure the average flux densities 
F v at 8, 24, and 70 jam wavelengths for the 441 regions. 
We do not consider the 160 /am band because its flux 
density relative to the 70 /am is much higher than is con- 
sistent with the DL07 models. We suspect that the 160 
/im flux is from cold dust that is not associated with 30 
Doradus but is in the sight line to the HII region. 

To ensure we are measuring the 8 and 24 /im flux 
densities only from dust and not starlight, we remove 
the starlight contribution at these wavelengths based on 
the 3.6/im flux density (which is almost entirely from 
starlight): 

F™(Sjim)=F v (%jim) - 0.232^(3.6//™) (4) 
F™(24/am)=F u (24/am) - 0.032F„(3.6/ira) (5) 

where the left-hand sides are the non-stellar flux at the 
respective wavelengths. The c oefficients 0.232 and 0.032 
are given in iHelou et~aT1 (|2004f ) . 

To account for the different spatial resolutions of the IR 
images, we convolved the 3.6, 8, and 24 /mi images with 
kernels to match the point-spread function of the 70 jam 
image. For this analysis, we employed t he convolution 
kernels and method described in §2.3 of Gor don et al.l 
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Fig. 7. — Plot of U versus PAH fraction qpAH- Arrows reflect 
upper-limits in the qpAH values, corresponding to the points out- 
side the grid in Fig. [6] 



(2008j). 

Figure [6] shows the resulting average IR flux ratios, 
(vF v )%l/(vF v ) 70 versus {vF u )% s /{i>F u )%l, of our regions. 
Overplotted are the DL07 model predictions for given 
values of #pah, the fraction of dust mass in PAHs, and 
U . Errors in our flux ratios are ~2.8% from a ~2% un- 
certainty in the Spitzer photometry. 

We interpolate the U-qpAU grid using Delaunay trian- 
gulation, a technique appropriate for a non-uniform grid, 
to find the U and #pah values for our regions. Figure 
plots the interpolated values of U versus #pah- Since the 
points with the smallest (uF u )g s / (vF^l values lie to the 
left of the U -#pah grid in Figure [6j we are only able to 
set upper limits of <7pah = 0.47% for them (marked with 
arrows in Fig.[7j). Thus, these regions produce the "wall" 
of points at #pah = 0.47% in the U versus #pah plot. 

We find that the PAH fraction spans roughly an order 
of magnitude, with values up to qpAU — 3.76%. U (and 
thus Uy) also varies significantly across 30 Doradus, with 
U w 91-7640, corresponding to u v w 6.6 x 10 -9 - 7.9 x 
10 -11 erg cm -3 . These ranges vary radially, with the 
largest U and smallest <?pah close to R136. One possible 
explanation for the <7pah radial dependence is that PAHs 
are destroyed more where the radiation field heating the 
dust is strong (e.g.. iGuhathakurta fc Draindl 19891). This 
result is consistent with the analyses of iPeeters et al.l 
(|2004f ). who showed that the ratio of PAH to far-IR (dust 
continuum) emission in Galactic HII regions is inversely 
correlated with the intensity of the UV field absorbed by 
the dust. 

We utilize the interpolated U values and Equation [3] 
to obtain the energy density u Vl and thus the pressure, 
of the dust-processed radiation field in the 441 regions. 

3.3. Warm Ionized Gas Pressure 

Next, we consider the pressure associated with both 
the warm HII gas and the hot X-ray gas. The warm ion- 
ized gas pressure is given by the ideal gas law, Phii = 
(n e + + ^He)fc^Hii? where n e , nn, and nn e are the 
electron, hydrogen, and helium number densities, respec- 
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Fig. 8. — Map of electron density n e of the warm gas across 30 
Doradus, calculated using the bremsstralung flux observed in the 
3-cm ATCA data. We find n e is fairly uniform, with a general 
range of n e ~ 100 — 500 cm -3 . 

tively, and Thii is temperature of the HII gas, which we 
assume to be the same for electrons and ions. If helium 
is singly ionized, then n e + nn + nn e ~ 2n e . The temper- 
ature of the HII gas in 30 Doradus is fairly homogeneous, 
with Thii = 10270 ± 140 K, based on the measurement of 
[O ill] ( A4959 + A5007) / A4363 ac ross 135 positions in the 
nebula (jKrabbe fc Co petti 2002); here, we adopt Thii = 
10 4 K. Since Thii is so uniform, the warm gas pressure 
is determined by the electron number density n e . We 
estimate n e from the average flux density F v of the free- 
free radio emission in each region (Eq. 5.14b, Rybicki & 
Lightman 1979): 



6.8 x 10 38 47rD 2 F u T^ 



1/2 



(6) 



where we have set the Gaunt factor ^ = 1.2. In the 
above relation, D is the distance to 30 Doradus (assumed 
to be D = 50 kpc) and V is the integrated volume of our 
regions. For V, we assume a radius of the HII region 
R=150 pc, and we calculate the volume by multiplying 
the area of our region squares by the path length through 
the sphere at the region's position. We measure F v of 
our regions in the 3.5-cm ATCA+Parkes image, since 
bremsstrahlung dominates at that wavelength. Figure [8] 
shows the resulting map of n e from these analyses. We 
find that the central few arcminute area of 30 Doradus 
has elevated electron densities, with values n e « 200—500 
cm -3 ; the location of these large electron densities corre- 
sponds to the two molecular clouds that form the " ridge" 
in the center of the nebula ([ Johansson et al.l [1998). The 
area outside the central n e enhancement has relatively 
uniform electron density, with n e « 100-200 cm -3 . In 
the Southwest of 30 Doradus where the supernova rem- 
nant N157B is located, we obtain elevated 3.5-cm flux 
densities, possibly because of a non-thermal contribution 
from that source. Therefore, the actual n e may be lower 
than the values we find in that region. 

Our warm gas electron densities are sim ilar to the val- 
ues obtained by llndebetouw et"aTl (|2009[ ) using Spitzer's 
Infrared Spectrograph. These authors used the ratio [S 
ill] A18.7 fim/[S ill] A33.4 fim to map n e across the cen- 
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tral ~2 7 of 30 Doradus. They also find enhancements in 
n e along the "ridge" . 

3.4. Hot Gas Pressure 

The hot X-ray gas arises from shock heating by stel- 
lar winds and supernovae, and the associated hot gas 
pressure is given by the relation Px = 1.9nxfcTx, where 
nx and Tx are the electron number density and temper- 
ature of the X-ray gas, respectively. The factor of 1.9 
arises from the assumption that He is fully ionized and 
that the He mass fraction is 0.3. As in our warm ion- 
ized gas calculation, we assume the electrons and ions 
have reached equipartition, so a single temperature de- 
scribes both. Since the hot gas can exist over a range 
of Tx, we measure both nx and Tx by modeling the 
bremsstrahlung spectrum at X-ray wavelengths. 

From the three Chandra observations, we extracted 
Chandra X-ray spectra from each region using the CIAO 
command specextract. Background spectra were also pro- 
duced from a circular region of radius «15" that is «2' 
East of 30 Doradus (the cyan circle in Figure [3]), and 
these were subtracted from the source spectra. Result- 
ing spectra were fit using XSPEC Version 12.4.0. Data 
were grouped such that a minimum of five counts were 
in each energy bin, and spectra from the three ACIS ob- 
servations of a given region were fit simultaneously to 
improve statistics (i.e., they were fit jointly, with more 
weight given to the data from the longer integrations). 
Around the edges of the HII region, the X-ray signal is 
weaker, so we combined adjacent regions to achieve suf- 
ficient counts for an accurate fit. 

Spectra were modeled as an absorbed hot diffuse gas 
in collisional ionization equilib rium using the XSPEC 
components phabs and mekal (|Mewe "eTatl [T985. 1986; 
Liedahl et al. 1995]). In this fit, we assumed a metallic- 
ity Z ~ Q.5Z(7), the value me asured in HII regions in the 
LMC (|Kurt fc Dufouri[T998h . For regions in the south- 
west of 30 Doradus with strong emission from the su- 
pernova remnant N157B, we added a power law compo- 
nent to account for the non-thermal emission from the 
SNR. We obtained good fits statistically, with reduced 
chi-squared values of 0.80-1.30 with 60-300 d.o.f. If a 
region's fit had reduced chi squared values outside this 
range or less than 60 d.o.f. we combined its spectra with 
those of adjacent regions to increase signal. The latter 
criterion was selected since we found generally that the 
shape of the bremsstrahlung continuum was not discern- 
able with less than 60 d.o.f. 

From our fits, we can estimate the electron number 
density nx of each region based on the emission mea- 
sure EM of our models. Emission measure is defined as 
EM = J n\dV. Thus, 



n x 



(7) 



where V is the integrated volume of our region (the same 
as used in Eq. [6]). Since we are interested in the contri- 
bution of the X-ray pressure to the global dynamics, we 
have divided the emission measure EM by the integrated 
volume of a region V in calculating nx, rather than the 
volume occupied by the hot gas. In the former case, the 
density nx goes as the filling factor / _1 ^ 2 ; in the latter 
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Fig. 9. — Map of X-ray gas temperature kTx (in keV) across 30 
Doradus. These values were obtained by modeling the Chandra 
X-ray spectra from each region. 
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Fig. 10. — Map of the hot gas electron density nx (in particles 
cm -3 ) across 30 Doradus. These values were obtained by modeling 
the Chandra X-ray spectra from each region, which output the 
best-fit emission measure EM. We converted EM to nx using 
Equation 7. 

scenario, nx oc / • / _1 ^ 2 = f 1 ^ 2 . If the filling factor of 
the hot gas is small, the thermal pressure of the bubbles 
may be high internally; however, the hot gas would be 
insignificant dynamically because it occupies a negligible 
volume and thus exerts little pressure on the material 
that bounds the HII region. Therefore, for our purposes 
of assessing the dynamical role of the hot gas, it is ap- 
propriate to use the integrated volue in calculating nx- 
Figure [9] shows a map of the best-fit temperatures kTx 
from the spectral modeling analyses. The X-ray gas tem- 
peratures are elevated in several areas, including in the 
Southwest (the bottom right of Fig. [9]), where the SNR 
N157B is located, and at the center near R136. Figure [T0l 
gives the map of the hot gas electron density across 30 
Doradus from our fits. We find a mean (nx) = 0.12 
cm -3 . The hot gas electron density is much less than that 
of the warm gas since many fewer electrons are heated 
to X-ray emitting temperatures (~ 10 T K) than to the 
moderate temperatures ~ 10 4 K of the warm gas. 
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Fig. 11. — All pressures versus radius from the center of R136. 
Regions with similar radii (denned as radii within 10% fractionally 
of each other) are binned to simplify the plot and make trends more 
readily apparent, and bars reflect the 1-a standard deviations in 
the pressures at the given radii. Generally, P^ r dominates at radii 
< 75 pc and follows a P^ir °c r ~ 2 relation (the blue solid line), 
whereas Phii dominates at larger distances from R136. Pjr and 
Px do not appear to contribute significantly. 

4. RESULTS 

Following the multiwavelength analyses and methodol- 
ogy of §3, we calculate the pressures associated with the 
direct stellar radiation pressure Pair, the dust-processed 
radiation pressure Pir, the warm ionized gas pressure 
Phii, and the hot X-ray gas pressure Px- Figure [TT] plots 
the results as a function of distance from the center of 
R136; data of similar radii (defined as radii within 10% 
fractionally of each other) are binned to simplify the plot 
and to make trends more readily apparent. By compar- 
ing the radial trends of the different components, we find 
that Pdir dominates at distances < 75 pc from R136, while 
Phii dominates at larger radii from R136. Additionally, 
Pir and Px do not appear to contribute significantly, al- 
though they are on the order of Pdi r at distances > 100 
pc from R136. 

As demonstrated by Figure HI we find that Phii > Px- 
The lack of pressure balance between these two compo- 
nents is consistent with our finding (see §5.1) that the 
X-ray gas does not remain adiabatic and trapped inside 
the shell. Instead, the hot gas is either leaking out or is 
mixing with cool gas and suffering rapid radiative losses 
as a result. In either case, the hot gas is likely to be flow- 
ing at a bulk speed comparable to its sound speed, and 
thus it will not have time to reach pressure equilibrium 
with the cooler gas that surrounds it before escaping the 
HII region. Alternatively, it may be that pressure bal- 
ance is established between the warm ionized gas and the 
ram pressure of the hot gas, whereas we have only mea- 
sured the thermal pressure. This picture is consistent 
with the anti-coincidence of th e warm and hot gas noted 
by pre vious X-ray work (e.g., Wang 1999; Townslev et al.l 
[20063). 

In Figure [I2j we give the maps of the four pressures 



across 30 Doradus for our 441 regions. Pdi r has a smooth 
profile due to its 1/r 2 dependence, while Phii is fairly 
uniform across 30 Doradus (as expected for a classical 
HII region). Compared to those components, Pir and Px 
have more variation throughout the source. Additionally, 
all the maps have significant enhancements in the central 
regions near R136; in the cases of Pir and Phii, the 
elevated pressures correspond to the molecular "ridge" 
in 30 Doradus (as seen in the CO contours in Figure [2]). 
Additionally, all except Pdi r have greater pressures in 
the regions near the SNR N157B (the bottom right of 
the maps). 

We can utilize the obtained pressures to estimate the 
total energy of each component. In particular, we mea- 
sure the total energy density u in a given radius bin of 
Figure [TT] and multiply by the volume of its shell (where 
we have set the shell thickness to the difference of the 
upper and lower bound radius of that bin). We con- 
vert pressures P to energy densities u using the rela- 
tions: P dir = w dir , Pir = ~i£ir, Phii = g^mi, and 
Px = |^x- Using this approach, we find the follow- 
ing total energies for each component: E^ Y = 5.1 x 10 53 
erg, E m = 1.7 x 10 53 erg, E mi = 2.8 x 10 53 erg, and 
Ex = 6.5 x 10 52 erg. Therefore, the direct and dust- 
processed radiation fields and the warm ionized gas con- 
tribute similarly to the energetics of the region, and every 
component is > 2 orders of magnitude above the typical 
kinetic energy of a single SN explosion. 

5. DISCUSSION 
5.1. Leakage of the Hot Gas 

As mentioned previously, the X-ray emission in 30 Do- 
radus arises from the shock-heating of gas to tempera- 
tures of ~ 10 T K by stellar winds and supernovae (SNe). 
These feedback processes eventually carve out large cav- 
ities, called bubbles and superbubbles, filled with diffuse 
X-ray emission. In Fig. [TTJ we demonstrated that the 
pressure associated with the hot gas Px is comparatively 
low relative to the other pressure components. Here, we 
explore the implications of this result in regard to the 
trapping/leakage of the hot gas. For this discussion, we 
will consider stellar winds only and ignore the contribu- 
tion by SNe; this assumption is reasonable given that 
the mechanical energy of one SN is on the order of the 
amount injected by winds over a single massive star's 
lifetime (Cas tor et al.lll975l ). This assumption is valid at 
the 0.5Z© of the LMC: simulations of a 5.5 x 10 4 M o star 
cluster in Starburst99 (jLeitherer et al.lfl999[ ) showed that 
the total wind luminosity decreased by roughly a factor 
of two from the solar to half-solar metallicity case. 

There are several competing theoretical models to 
account for the X-ray lumin osity in bubb l es an d su- 
perbubbles. The m odels of iCastor et al.l (|1975l ) and 
IWeaver et all (|1977[ > assumes that the shock-heated gas 
is completely confined by a cool shell expanding into a 
uniform- density ISM. An alternative theory proposed by 
IChevalier fc Clegel fl985) ignores the surrounding ISM 
and employs a steady-state, fre e-flowing wind. Recently, 
lHarper-Clark fc Murray! (2009) introduced an intermedi- 
ate model between these two, whereby the ambient ISM 
is non-uniform. In this case, only some of the hot gas 
can escape freely through the holes in the shell. 

The fraction of hot gas confined by the shell directly 
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Fig. 12. — Maps o f th e four pressure components across 30 Doradus. All four are on the same color scale to enable visual comparison. 
Consistent with Fig. 1111 Pdi r dominates in the central few arcminutes, while the Phii dominates at larger distances from R136. 



determines the hot gas pressure on the shell as well as the 
X-ray luminosity within the bubble. If the shell is very 
porous, the shock- heated gas will escape easily, the wind 
energy will be lost from the bubble, and the associated 
pressure and luminosity will be low. By comparison, a 
more uniform shell will trap the hot gas, retain the wind 
energy within the bubble, and the corresponding X-ray 
pressure and luminosity will be much greater. As such, in 
the latter case, the shocked winds could have a significant 
role in the dynamics of the HII region. We note that the 
warm gas is not able to leak similarly because its sound 
speed is less than the velocities of the shells (20-200 km 
s- 1 : IChu fc KennlcTrttlflQQl . 

To assess whether the hot gas is trapped inside the 
shell and is dynamically important, we measure the ratio 
of the hot gas pressure to the direct radiation pressure, 
/tra P ,x = Px/Pdir, and compare it to what / tra p,x would 
be if all the wind energy was confined. We can calculate 
the t rapped-wind value us i ng; the wind- luminosi ty rela- 
tion (jKudritzki et al.l 119991 : iRepolust et al.ll2004f ). which 
indicates that the momentum flux carried by winds from 
a star cluster is about half that carried by the radia- 
tion field if the cluster samples the entire IMF. Written 
quantitatively, 0.5Lt> o i/c = M w v w , where M w is the mass 
flux from the winds that launched at a velocity v w . The 
mechanical energy loss L w of the winds is then given by 



L w -M w v w 



T 2 



8M w c 



2' 



(8) 



and the mechanical energy of the winds is simply E w = 
L w t, where t is the time since the winds were launched. 
Putting these relations together, the trapped X-ray gas 
pressure Px,t is 



Px; 



2E W 
3Vkn 



Lit 



1QttM w c 2 RI u 



(9) 



where Vmi is the volume of the HII region. 

Given that P dir = L b oi/(47ri% n c), then / tra p,x is 



«/trap,X 



^bol 



(10) 



4M w cR m i 4M w cv sh 

where we have set Rmi/t = ^sh, the velocity of the ex- 
panding shell. Finally, we put M w in terms of Lboi and 



u w , so that Eq. fTUI reduces to 

«/trap,X = 



2^ s h 



(11) 



We use the above equation to obtain an order-of- 
magnitude estimate of /t ra p,x if all the wind energy 
is confined by the shell. We assume a wind velocity 
v w ~ 1000 km s _1 (the escape velocity from a 06 V star; 
a reasonable order-of-magnitude estimate, since 03 stel- 
lar winds are faster and WR winds would be slower than 
this value). If we set ~ 25 km s _1 (the expansion 
velocity over 30 Doradu s given by optical spectroscopy; 
IChu fc KennTc7rttl[l99l . then / trap ,x ~ 20. 

We can compare this /t ra p,x to our observed values for 
the regions closest to the shell (the ones along the rim 
of our 441 squares in Fig. [5j); Figure [13] shows the his- 
togram of our observed /t ra p,x values. We find a mean 
and median /t ra p,x of 0.30 and 0.27, respectively, for our 
outermost regions. Over 30 Doradus, the highest val- 
ues of /trap,x are near the supernova remnant N157B 
in the southwest corner of 30 Doradus (see Figure [H]), 
where hot gas is being generated and has not had time 
to vent. Other locations where /t ra p,x is elevated are re- 
gions with strong X-ray emission and weak Ha emission. 
Morphologically, these areas could be where the hot gas 
is blowing out the 30 Doradus shell. 

The observed /t ra p,x values are 1-2 orders of magni- 
tude below what they would be if the wind was fully con- 
fined. As a consequence, we find that Px of our regions is 
too low to be completely trapped in the HII region (the 
Castor et al. model), and the X-ray gas must be leaking 
through pores in the shell. This result is consistent with 
the Harper- Clark & Murray model of partial confinement 
of the hot gas, and the weakness of Px relative to Pdi r 
suggests the hot gas does not play a significant role in 
the dynamics of the HII region. 

We note here that our rim regions in this analysis are 
~ 70-130 pc from R136, which is less than the estimated 
radius of Rmi = 110 — 185 pc. Therefore, our /t rap ,x 
values are lower limits, and the true /t rap ,x at the shell 
may be greater by a factor of a few. Nonetheless, the 
conclusions would remain the same. 

An alternative explanation for the weak X-ray lumi- 
nosity is that the hot gas mixes with the cool gas, and 
the hot gas temperature is lowered enough so that the 
gas can cool efficiently. In that case, the energy is still 
lost from the system, via radiative cooling instead of the 
escape of the X-ray emitting material. Far ultraviolet 
spectroscopy is necessary to determine the level of mix- 
ing between the gas components. 
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Fig. 13. — Histogram of /trap,x = Px/Pdin the ratio of hot gas 
pressure to the direct radiation pressure, for the regions that are 
along the rim of our 441 squares in Fig. [5] We find that the mean 
/trap,x is 0-30 and the median is 0.27, far below the values expected 
if the hot gas is completely confined in 30 Doradus (/trap,x ~ 20; 
see text). This result is evidence that the hot gas is leaking out of 
the shell. 




Fig. 14. — Ha (green) and soft X-ray (blue) images with the 
region grid overplotted, where the regions are color-coded based 
on their / tr a P ,x = Px/Pdir values: dark blue = / tr ap,x < 0.2; 
cyan = 0.2 < / tra p,x < 0-3; yellow = 0.3 < / tra p,x < 0-5; red = 
/tra P ,x > 0.5. The highest values of /t ra p,x occur near the SNR 
N157B as well as other areas with strong X-ray emission and weak 
Ha emission. 

In this paper, we have defined the radiation pressure 
as related to the energy and momentum flux of the light 
radiated by the stars in 30 Doradus. Alternatively, ra- 
diation pressure could be characterized as the force per 
unit area exerted by the radiation on matter. The two 
cases produce divergent results regarding the radial de- 



pendence of Pdir- In particular, the former case has large 
Pdir close to the star cluster and a decline in Pdir with dis- 
tance from the center. By contrast, the latter predicts 
small Pdir in the HH-region interior, and Pdi r becomes 
significant near the neutral shell where the radiation will 
be absorbed (see Appendix §B). 

Each definition of radiation pressure reveals distinct 
information about an HII region. When considering the 
global dynamics of expansion of an HII region, it is nec- 
essary to characterize Pdir as the energy density of the 
radiation field, since that definition reflects the total en- 
ergy and momentum budget available to drive motion. 
Alternatively, measurement of the force exerted by ra- 
diation on matter facilitates a local estimate of the in- 
ternal density distribution of an HII region. As we are 
interested principally on the dynamical role of radiation 
pressure, we have adopted the former definition of Pdi r 
in this paper. 

5.3. HII Region Dynamics 

In we found that the direct radiation pressure 
Pdir dominates over the ionized gas pressure Phii at 
radii < 75 pc, implying that the radiation has played 
a role in the dynamics in 30 Doradus. Significant ra- 
diation pressure alters th e properties of an HII region 
(e.g., the density profile: iDraind [201 Of ) and causes the 
expansion to proceed differently than in a classical HII 
regio n with ionized gas-dri v en ex pansion. In partic- 
ular, iKrumholz fc Matznerl ([2009) demonstrated that 
radiation-driven expansion imparts more momentum to 
the shell, accelerating the expansion at early times rel- 
ative to that of gas-driven expansion. Indeed, an ad- 
ditional force must have dominated at early times in 
30 Doradus since the s hell velocity v s h ~ 25 km s _1 
(jChu fc Kennicuttl Il994f ) is too fast to have been gas- 
driven alone because the HII gas sound speed is c s ~ 10 
km s _1 . 

To determine the characteristic radius r c h where the 
HII region shell transitioned from radiation-pressure 
driven to gas-pressure driven, we can set these pressure 
terms at the shell equal and solve for r c h- Broadly, the 
pressures at the shell have different dependences with 

rj 3/2 

the shell radius rmi- Pdir oc r mi and Phii cx r m { . 
Setting them equal and so lving for r c h (Equation 4 in 
IKrumholz fc Matznerl I2009L the "embedded case" ) , we 
find 



Tch ~ 12^\2.2k B T mi J tmp,t 



(12) 



where a B is the case-B recombination coefficient and 
eo = 13.6 eV, the photon energy necessary to ionize hy- 
drogen. The dimensionless quantity <\> accounts for dust 
absorption of ionizing photons and for free electrons from 
elements besides hydrogen; <j> = 0.73 if He is singly ion- 
ized and 27% of photons are absorbe d by dust (typical for 
a gas -pressure dominated HII region; McKee fc Williams! 
119971 ). The /trap, tot represents the factor by which ra- 
diation pressure is enhanced by trapping energy in the 
shell through several mechanisms; the trapped hot gas 
/trap,x calculated in §5.1 is one co mpon ent that can con- 
tribute to /trap, tot (as discussed in j35.1D. Here, we adopt 
/trap, tot = 2, as in IKrumholz fc Matznerl (j2009[ ). Lastly, 
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ip is the ratio of bolometric power to the ionizing power 
in a cluster; we set ijj = 3.2 using the (S) / (M*) and the 
(L)/(M*) relations of lMurrav fc Rahmanl (|201Q[ ). 

Putting all these terms together, we find r c h ~ 33 pc. 
Physically, this result means that early in the expansion 
before it reached a radius of 33 pc, 30 Dor's dynamics 
could have been radiation-pressure dominated, and it has 
since become gas-pressure dominated. Alternatively, it 
is possible that the hot gas pressure dominated at early 
times and has become weaker as the HII region expands. 

The radiation-driven or hot gas driven expansion at 
early times in 30 Doradus would have facilitated the 
expulsion of gas from the central star cluster. In par- 
ticular, since the warm gas sound speed (~10 km s _1 ) 
is less than the escape velocity of R136 (~20 km s _1 , 
given a mass \ M = 5.5 x 1O 4 M in a radius R = 1 pc; 
iHunter et al.|[l995[ ). an alternative mechanism is neces- 
sary to remove the gas and r egulate star form ation (e.g., 
iKrumholz fc Matznerl l200l iFall et al.ll20lol ). We con- 
clude that the radiation pressure or hot gas pressure 
likely played this role in 30 Doradus, decreasing the avail- 
able mass to make new stars and slowing star formation 
in the region. 



6. SUMMARY 

In this paper, we have utilized mult i- wavelength (ra- 
dio, infrared, optical/UV, and X-ray) imaging to assess 
the role of several stellar feedback mechanisms in the gi- 
ant HII region 30 Doradus in the LMC. In particular, 
we have measured observationally the pressures associ- 
ated with possible sources of energy and momentum to 
drive the dynamics of the region: the direct radiation 
from stars, the dust-processed infrared radiation field, 
the warm ionized gas from massive stars, and the hot 
gas shock- heated by stellar winds and supernovae. We 
have exploited the high-resolution images of 30 Doradus 
to map these pressure components in 441 square regions, 
with dimensions of 35" x 35". We have found that the di- 
rect radiation pressure from stars dominates at distances 
less than 75 pc from the central star cluster, whereas 
the warm ionized gas pressure dominates at larger radii. 
By contrast, the hot gas pressure and the dust-processed 
radiation pressure do not contribute significantly, indi- 



cating these components are not dynamically important. 
However, we cannot rule out that the hot gas pressure 
dominated at early times and has become weaker with 
the HII region expansion. 

We have discussed two implications of our results: the 
partial confinement of the hot gas and the dynamical role 
of radiation pressure in 30 Doradus. First, the weakness 
of the X-ray gas pressure relative to the direct radiation 
pressure suggests the hot gas is only partially confined 
and is leaking out of the pores in the HII shell. Secondly, 
the significant radiation pressure near the star cluster in- 
dicates that radiation pressure may have driven the ex- 
pansion of the HII shell at early times. This result sug- 
gests observationally that radiation pressure may be dy- 
namically important in massive star clusters, reinforcing 
that radiation pressure is a viable mechanism to remove 
gas from HII regions and to regulate star formation. In- 
deed, if NGC 2070 was more massive, the radiation pres- 
sure could even expel gas at high enough v elocities to 
launch a galactic wind (jMurrav et al.ll2010al ). 

The work presented here is a first step to measure ob- 
servationally the relative role of stellar feedback mecha- 
nisms in star- forming regions. Although we have applied 
our techniques to one source, 30 Doradus, our methods 
to extract dynamical information from multi- wavelength 
images can be applied to other sources as well. Con- 
sequently, we plan to perform these analyses on all the 
HII regions in the LMC with available data to develop a 
broad observational understanding of these stellar feed- 
back mechanisms and their role in regulating star forma- 
tion. 
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APPENDIX 
DEPROJECTING THE STARS IN 30 DOR 

Since the stars are viewed in projection, the actual distance r to a star from the R136 center is observed as a 
projected distance i/j (see Figure [15]). Therefore, we calculate the direct radiation pressure for two scenarios: one case 
assuming the stars lie in the same plane (i.e., r = ip; Pdir) and another case where we attempt to "deproject" the stars 
positions (i.e., r = vV 2 + z<2 i where z is the line-of-sight distance to the star from the sphere's midplane; Pdir, 3D )• The 
direct observable is the projected surface brightness fi (in units of erg cm -2 s _1 ) in an annulus, and it is a function of 
ip. The luminosity density (in units of erg cm -3 s _1 ) C{r) = C{y/^ 2 + z 2 ) is then related to /i(V0 by 



r ^R 2 -*p 2 

/i(V0 = 2 / C(y^ 2 ~T^)dz. 
Jo 



If we put this integral in terms of r, we obtain the relation 

tip JR 



rR ftp 

H(tf) = 2 r(r 2 - i/j 2 )- 1/2 £(r)dr = / K(r^)C(r)dr 

Jib JR 



(Al) 



(A2) 
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Fig. 15. — Diagram explaining how projection effects may lead to erroneous measurement of the actual distance r of a star to the star 
cluster center, i. Face-on view of an HII region of radius R. The projected distance from the star to the center is tp. ii. View from above 
of the same HII region. In this case, it is apparent that the star does not lie in the midplane of the sphere, and the actual distance is 
r = \f^ 2 + z 2 , where z is the line-of-sight distance to the star from the sphere's midplane. 
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Fig. 16. — Comparison of Pdir, 3D an d -Pdir- Top: The two pressure components versus distance R from R136. The two have similar radial 
dependence. Bottom: The fractional difference (Pdir — Pdir,3D)/Pdir for all the points in the top panel. At radii > 20 pc from R136, the 
fractional difference is small (0.1%-3.0%), while at radii < 20 pc, the fractional difference is greater (10%-60%). 



where K(r, ijA = — 2r(r 2 - ^ 2 ) -1 / 2 . 

Equation IA2I is a Volterra equation of the first kind, and we solve for C(r) explicitly for annuli beginning at some 
radius R and taking N uniform steps of size h inward to r\ — R — hN. In this case, we selected a radius R = 200 
pc « 825.5 // to ensure the entire nebular volume was included. Additionally, we chose a step size of h = l"and went 
inward to T\ = 17.5"(so N = 808 steps), the radius of the HST PC image. 

Figure [16] (top panel) plots the resulting Pdir, 3D (and Pdi r for comparison) versus the distance R from R136. The 
bottom panel gives the fractional difference (Pdir — Pdir,3D)/Pdir for all the points in the top panel. The fractional 
difference between Pdir, 3D and Pdi r is small (~0.1%-3.0%) for regions ^20 pc from R136, and becomes larger (~10%- 
60% for radii < 20 pc. Despite these greater fractional differences at smaller radii, Pdir, 3D would still dominate over 
the other pressure components in Fig. [11] at distances < 75 pc from R136. Additionally, the small fractional differences 
at distances > 20 pc confirms that our values of Pdi r near the shell are accurate. 

We note in the above calculation, we necessarily assumed that the luminosity in a shell is spherically symmetric. 
The surface brightness in 30 Doradus is not actually symmetric though (since stars and star clusters are distributed 
non-uniformly around the nebula), so our estimated Pdir, 3D values are an approximation of the true, deprojected 
radiation pressure. Nonetheless, the small differences between Pdir, 3D and Pdi r indicate that uncertainty in the star 
position along the line of sight does not qualitatively affect our results. 
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r/ R 

Fig. 17. — Plot of /, the force per unit volume on matter from radiation, versus fractional radius r/R of the HII region. The force density 
/ peaks near the HII shell and is several orders of magnitude less interior to the shell. 



ON THE DEFINITION OF RADIATION PRESSURE 

An alternate definition of radiation pressure than the one used in this paper is to characterize Pdi r as the force per 
unit area exerted by the radiation on matter. This case predicts small Pdir in the HH-region interior, where the density 
is small, and Pdir only becomes significant near the neutral shell where the radiation is absorbed. To demonstrate 
this effect, we calculate /, the force per unit volume on matter from radiation, as a funct ion of radius in an idealized 
HII region using Version 08.00 of the photoionization code cloudy (jFerland et al.|[l998l ). Assuming photoionization 
balance, we have 

S(r) 

a B n e n p = -^^a m n m , (Bl) 

where a B = 2.6 x 10 -13 cm 3 s _1 is the case-B recombination coefficient of hydrogen at 10 4 K, n e , n p , nni are the 
electron, proton, and HI number densities, S is the ionizing photon luminosity (in photons s _1 ) passing through a 
shell at a given radius r, and ctri is the H ionization cross-section. 
The force density / is given by 

npF ^dustP^ . ftHipP /d m 
J = 1 , {HI) 

c c c 

where n is the opacity per unit mass (from dust ftdust or fr° m neutral hydrogen ^m), P is the local mass density, and 
F is the total flux. Given ftdustP = ^H^dust (where n B = ^hi + ^hii) and ftmP = ^hi^hi, 

S(r){hv) a B n e n p (hv) ( n n a dnst \ 

/ — ; — o — (^HCTdust + nmcrm) = 1 H • (B3) 

47rr z c c \ nni <thi / 

The first term represents the force of ionizing photon absorption by H atoms and the second term is the force of ionizing 
photons on the dust. Here, (his) is the mean energy of the ionizing photons, assuming the force in the radiation field is 
from ionizing photons only. In the following calculation, we set (his) = 15 eV, a value typical of an O star. Similarly, 
we assu me the force on dust is from ionizing photons, and we adopt a dust cross section at 15 eV, <7d U st = 1-0 x 10 -21 
cm 2 /H (Weingar tner fc Draindl200lh . In the above expression, we set am = 6.3 x 10 -18 cm -2 . 

The radial dependence of / comes from the density profiles, n e (r) and n p (r). To estimate these parameters as 
a function of radius in CLOUDY, we utilize the OSTAR2002 stellar atmosphere model for a metallicity Z/Zq = 0.5 
(jLanz fc Hubenvl l2003) and stellar temperature of T* = 35000 K. In our analysis, we set our idealized HII region to 
have S = 6.9 x 10 49 photons s _1 and nn = 10 cm -3 . Additionally, we include the effects of interstellar grains in the 
calculation. 

Figure [T71 plots the resulting / as a function of fractional radius in the HII region. The force per unit volume of the 
radiation increases many orders of magnitude from the center to the edge of the HII region, with a drastic jump in / 
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at r/R> 0.9. This plot contrasts Figure [TTj where Pdir falls off as 1/r 2 . 
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